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Abstract 


We  report  on  the  study  of  Josephson  junction  two-cell  flux  qubits  and  phase  qubits. 
We  have  made  progress  with  fabrication  of  two-cell  flux  qubits  based  on  small  Al-A10x- 
A1  Josephson  tunnel  junctions.  We  succeeded  in  preparing  junctions  as  small  as  0.09  x 
0.24  /im2.  We  describe  newly  developed  sample  design  which  implements  a  “silent”  two¬ 
cell  qubit  containing  a  fluxon  trapping  loop  as  a  phase  shifter.  This  design  was  submitted 
for  fabrication  at  MIT  Lincoln  Lab  foundry.  We  characterized  several  chips  fabricated 
according  to  this  design  at  MIT  Lincoln  Laboratory  and  observed  a  switching  of  the  qubit 
from  one  classical  state  to  the  other.  A  new  readout  scheme  which  uses  fast  current  pulses 
is  under  development. 

The  states  of  superconducting  quantum  bits,  as  the  Josephson  flux  and  phase  qubits, 
are  distinguished  by  magnetic  flux  generated  by  a  persistent  current  in  a  superconducting 
loop.  This  flux  is  read  out  by  measuring  the  switching  current  of  an  inductively  coupled 
SQUID  detector.  As  the  readout  circuit  contributes  decoherence,  its  coupling  to  the  qubit 
should  be  made  as  weak  as  possible.  In  the  second  part  of  this  report,  we  demonstrate  a 
new  data  evaluation  method  which  maximizes  the  measurement  contrast  of  the  two  qubit 
states,  thus  allowing  to  significantly  decrease  the  coupling.  The  method  is  tested  experi¬ 
mentally  by  measuring  high-fidelity  Rabi  oscillations  in  phase  qubits  operated  beyond  the 
single-shot  resolution  limit.  This  technique  is  promising  for  multi-qubit  circuits  since  it 
allows  to  read  out  the  states  of  several  qubits  by  a  single  weakly  coupled  SQUID. 

Our  current  work  on  Josephson  phase  qubits  is  focussed  on  the  study  of  temperature 
effects  on  decoherence.  Using  samples  fabricated  from  different  materials  and  origins,  we 
measured  Rabi  and  Ramsey  oscillations  in  a  wide  range  of  temperatures.  Unexpectedly, 
the  decoherence  times  extracted  from  the  data  show  only  a  very  weak  dependence  on 
temperature.  We  find  that  coherent  operation  of  the  qubits  remains  possible  as  long  as  the 
thermal  energy  fee  T  is  smaller  than  the  energy  difference  between  the  qubit  eigenstates. 
Since  the  latter  is  determined  by  the  circuit  design,  it  makes  possible  operating  currently 
available  superconducting  qubits  at  temperatures  of  300  mK  and  even  higher. 


1  Two-cell  flux  qubits 

The  flux  qubits  use  junctions  enclosed  in  superconducting  loops  as  quantum  coherent 
systems.  Single-cell  flux  qubits  have  been  previously  studied  at  MIT  and  Delft.  These 
devices  have  already  shown  quantum  superposition  of  macroscopic  persistent  current  states 
and  quantum-coherent  oscillations  in  direct  time-resolved  measurements.  We  are  focusing  on 
two-cell  flux  qubits  consisting  of  four  identical  Josephson  tunnel  junctions.  These  devices 
have  an  advantage  of  having  a  control  over  the  barrier  hight  between  two  potential  wells. 

1.1  “Silent”  two-cell  qubit 

We  developed  a  new  scheme  for  two-cell  flux  qubits  which  we  call  a  “silent”  qubit.  Its 
main  idea  can  be  seen  from  Fig.  [l]  Having  a  flux  quantum  trapped  permanently  in  the  inner 
loop,  the  qubit  has  the  advantage  of  using  smaller  range  of  magnetic  fields  for  the  symmetric 
flux  control.  The  trapped  flux  serves  as  7r-shifter  of  the  phase  and  contributes  equally  to  the 
phase  change  in  both  loops.  Thus  one  does  not  need  to  apply  large  external  flux  and  the  qubit 
is  operated  close  to  zero  symmetric  field.  In  order  to  avoid  steep  dependence  of  qubit  level 
splitting  on  magnetic  field,  we  can  make  junctions  in  outer  qubit  loop  smaller  as  suggested 
recently  by  Dr.  Stanford  Yukon. 

A  two  cell  flux  qubit  contains  four  Josephson  junctions  of  two  different  sizes.  Two  outer 
Josephson  junctions  numbered  #3  and  #4  have  smaller  area  than  two  inner  ones  numbered 
#1  and  #2.  We  denote  by  p  the  ratio  between  their  critical  currents  and  thus  their  Josephson 


2 


A 

k 

9i) 

ct 

<P3  _  | 
xt  cy 

\pa  "T 

-6 

— x — 

92) 

ct 

Figure  1:  A  “silent”  two-cell  flux  qubit 
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The  potential  energy  U  of  the  qubit  is  periodic  as  a  function  of  the  average  phase  x  = 
0.5(cpi  +  P2)  on  junctions  #1  and  #2: 


U  =  —  2  cos  x  +  t  cos(2y  —  4>a)  ■  (2) 

Here  r  =  — 2/ocos(0s),  while  <f>8  =  and  4>a  =  ^L~^R  are  the  symmetric  and  anti¬ 

symmetric  parts  of  the  flux  in  the  two  cells,  respectively.  Each  period  of  U(x)  has  two 
minima.  For  p  =  r/ 2  we  can  operate  this  qubit  at  half  frustration  in  each  cell.  With  a  single 
flux  quantum  d>o  frozen  in  the  inner  loop  this  would  require  zero  or  (in  non-ideal  case)  a  very 
small  external  magnetic  field.  The  inner  loop  thus  acts  as  passive  7r-shifter.  The  potential 
([2])  has  two  minima  associated  with  the  circulating  currents  in  the  outer  loops,  as  shown  in 
Fig.  m  At  the  optimal  working  point  of  r  =  0.76,  the  current  through  the  central  line  changes 
from  state  (a)  to  state  (b)  by  approximately  2/c3’4^  =  2 pl^'2\ 


(a) 


(b) 


Figure  2:  Two  states  of  the  two-cell  flux  qubit  shown  in  Fig.  [l]  The  arrows  show  the  flow 
directions  of  the  persistent  currents. 

In  quantum  regime,  at  the  optimal  operating  point  the  mean  current  through  every  Joseph- 
son  junction  is  zero.  Taking  into  account  that  we  work  near  zero  external  magnetic  field,  it 
means  that  at  the  operating  point  the  “silent”  qubit  is  practically  decoupled  from  the  fluctu¬ 
ations  caused  by  the  flux  control  and  readout  lines. 

1.2  Nb/AlO^/Nb  junction  based  qubits 

1.2.1  Samples  fabricated  at  MIT  Lincoln  Laboratory 

At  the  end  of  December  2006  we  submitted  a  new  design  to  MIT  Lincoln  Laboratory 
foundry.  We  designed  a  readout  circuit  to  measure  the  “silent”  qubit  described  above.  This 
circuit  contains  an  asymmetric  three-junction  dc-SQUID  inductively  coupled  to  the  qubit. 
The  SQUID  has  sufficient  flux  sensitivity  at  zero  external  magnetic  field.  The  qubit  is  placed 
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inside  the  readout  SQUID  loop,  which  is  coupled  asymmetrically  to  the  two  cells  of  the  qubit. 
The  flux  through  the  SQUID  loop  changes  when  the  qubit  switches  between  two  persistent 
current  states.  The  targeted  critical  current  Icq  of  inner  qubit  junctions  is  about  0.5  //A  and 
capacitance  Cq  ~  5.8  fF.  The  outer  qubit  junctions  are  designed  2.5  times  smaller  then  the 
inner  ones.  The  calculated  level  splitting  of  the  qubit  is  around  oj oi  ~  1.5  GHz  at  the  working 
point  r  =  0.8.  We  optimized  self  and  mutual  inductances  of  the  qubit  cells  and  the  readout 
SQUID  to  have  the  change  of  the  flux  through  the  SQUID  loop  induced  by  qubit  switching 
between  two  states  to  be  about  AT  ~  0.002  To-  Sizes  of  junctions  were  optimized  for  three 
announced  critical  current  densities  of  500  A/cm2,  1  kA/cm2  and  2  kA/cm2.  A  fragment  of 
the  submitted  layout  is  presented  in  Fig.  [3j  In  this  circuits  we  intend  to  use  biasing  resistors 
on  chip  (shown  in  green).  We  have  placed  on  chip  three  different  types  of  control  lines:  (i)  an 
antisymmetric  line  to  change  antisymmetric  part  of  the  flux  through  the  qubit  cells  (in  Fig.  [3] 
it  is  situated  to  the  right  of  the  qubit),  (ii)  a  symmetric  line  to  control  the  working  point  r 
(to  the  left  of  the  qubit  in  Fig.  [3]),  and  (iii)  two  lines  for  biasing  the  readout  SQUID  (below 
the  qubit  in  Fig.  [3]) . 


Figure  3:  A  part  of  the  design  submitted  to  the  MIT  Lincoln  Laboratory. 

In  October  2007  we  received  two  batches  of  Nb/AlO^/Nb  junction  chips  fabricated  ac¬ 
cording  to  our  design  at  MIT  Lincoln  Laboratory.  Every  chip  contains  three  sets  of  structures 
designed  for  critical  current  densities  of  500  A/cm2,  lkA/cm2  and  2kA/cm2.  Every  set  con¬ 
sists  of  three  qubit  circuits  with  different  coupling  strength  between  the  readout  SQUID  and 
the  qubit  and  four  single  Josephson  junctions  which  have  the  same  sizes  as  junctions  in  the 
qubit  and  the  readout  SQUID.  We  use  them  for  testing  purposes.  In  qubit  circuits  we  use 
junctions  of  three  different  areas.  To  check  the  process  reproducibility  at  a  small  scale,  two 
test  junctions  have  the  size  equal  to  the  size  of  the  smallest  qubit  junctions. 

Figure  [4]  shows  an  SEM  micrograph  of  a  two-cell  flux  qubit  together  with  its  readout 
SQUID. 

The  fabrication  process  is  characterized  by  a  process  bias,  which  determines  a  reduction 
of  the  electrical  junction  sizes  form  the  drawn  sizes.  The  average  process  bias  in  the  batch 
delivered  by  Lincoln  Laboratory  is  about  0.25  gm.  This  value  was  taken  into  account  when 
designing  the  qubits.  Larger  process  biases  result  in  increasing  the  maximal  barrier  hight 
and  decreasing  of  the  qubit  flux  signal.  Smaller  process  biases  result  in  decreasing  of  the 
maximal  barrier  hight,  this  makes  qubit  operating  difficult.  However,  there  is  a  variation 
across  the  wafer  and  from  wafer  to  wafer.  For  the  wafer  DSM3_07_3_8  the  process  bias  is  to 
0.23  ±  0.03  fi m.  The  dies  with  the  process  bias  close  to  0.25  /urn  had  smaller  yield  of  junctions. 
The  yield  for  those  chips  is  seven  working  junctions  out  of  eleven  drawn,  the  median  yield  over 
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O  Josephson  junctions 


Figure  4:  An  SEM  micrograph  of  a  sample  fabricated  at  MIT  Lincoln  Laboratory. 

the  wafer  is  nine  working  junctions.  The  yield  estimate  includes  the  drawn  0.2  /jm  and  0.3  /mi 
Josephson  junctions.  The  smallest  junction  reliably  yielded  has  the  drawn  size  of  0.4  /mi, 
while  the  smallest  qubit  junction  has  the  drawn  size  of  0.45  /mi.  For  the  wafer  DSM4_07_11_4 
the  process  bias  is  0.21  ±  0.04  /mi  and  the  median  yield  is  eight  working  junctions  out  of 
nine  drawn.  The  smallest  junction  has  the  drawn  size  of  0.4  /mi.  The  chip  B4  is  the  only  die 
with  the  process  bias  equal  to  0.25  /mi  on  this  wafer.  There  are  several  dies  with  process  bias 
0.26  /mi  and  0.27  /mi  on  the  wafer,  but  the  expected  qubit  signals  on  these  chips  might  be  too 
week,  as  the  critical  current  density  on  these  dies  is  around  440  A/cm2.  Thus  only  few  chips 
from  the  available  batch  were  expected  to  be  close  to  the  targeted  parameters. 


1.2.2  Measurements 

The  test  Josephson  junctions  showed  high-quality  current-voltage  characteristics,  which 
are  presented  in  Fig.  [5j  By  scaling  these  curves  we  can  calculate  sizes  of  the  junctions.  The 
small  qubit  junctions  have  area  of  about  0.03 /mi2,  the  large  qubit  junctions  are  of  about 
0.08 /mi2  and  the  large  SQUID  junctions  have  the  area  of  about  0.16 /mi2.  These  values 
are  close  to  designed  values,  which  were  0.03 /mi2,  0.07/mi2  and  0.14 /mi2  respectively,  and 
suitable  for  qubit  operating.  One  can  readily  see  that  the  spread  in  critical  currents  of  two 
smallest  qubit  junctions  is  negligible. 

Figure  [6]  shows  measured  current- volt  age  characteristics  of  readout  SQUIDs  on  the  chip 
C4  from  the  wafer  DSM3_07_3_8  and  on  the  cpip  B4  from  the  wafer  DSM4_07_11_4.  Despite 
the  high  quality  of  the  test  junctions,  the  critical  current  of  the  readout  SQUID  is  much 
below  the  expected  value.  There  is  no  peculiarity  around  3  mV,  which  should  correspond 
to  2 A/e  ,  where  A  is  the  superconducting  energy  gap  of  niobium.  We  find  a  step  around 
lmV  for  the  wafer  DSM3_07_3_8  and  0.4  mV  for  the  wafer  DSM4_07_11_4,  which  can  most 
probably  be  due  to  difference  of  energy  gaps  in  top  and  bottom  niobium  electrodes.  We 
measured  six  chips  from  both  wafers,  and  all  15  measured  readout  SQUIDs  had  the  current- 
voltage  characteristics  similar  to  ones  presented  in  Fig.  [6}  A  couple  of  SQUIDs  showed  clear 
indications  of  shorts,  but  they  might  have  been  destroyed  by  an  electro-static  discharge  during 
mounting  procedure. 
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Figure  5:  Current-voltage  characteristics  of  the  test  junctions  on  chip  C4  from  the  wafer 
DSM3_07_3_8  at  300  rnK. 

We  observed  a  modulation  of  the  SQUID  critical  current  with  applied  magnetic  flux.  The 
dependence  of  the  SQUID  critical  current  on  symmetric  control  flux  induced  by  control  line 
Is  for  the  readout  SQUID  on  the  chip  B4  from  the  wafer  DSM4_07_11_4  is  presented  in  Fig.  [Tj 
The  plot  shows  that  the  SQUID  has  non-zero  sensitivity  at  zero  external  magnetic  field.  This 
allows  us  to  read  the  qubit  state  out  around  zero  symmetric  control  flux. 

We  believe  that  the  zero-field  response  of  the  circuit  is  due  to  an  odd  number  of  flux 
quanta  trapped  in  the  passive  7r-shifter  loop  while  cooling  down.  We  made  three  attempts 
to  trap  flux  inside  the  passive  7r-shifter  loop.  Making  further  tries  turned  impossible  due  to 
rapid  degradation  of  the  readout  circuit.  We  saw  a  shift  of  the  dependence  of  the  SQUID 
critical  current  on  the  control  current  Is  after  every  attempt,  which  means  that  we  trapped 
a  different  number  of  quanta.  However,  every  time  we  trapped  an  odd  number  of  quanta,  as 
we  saw  the  qubit  switching  at  zero  applied  symmetric  flux. 

Figure  [8]  shows  the  dependence  of  the  readout  SQUID  on  chip  B4  from  the  wafer 
DSM4_07_11_4  critical  current  on  antisymmetric  control  flux  at  zero  applied  symmetric  flux. 
The  dependence  is  non-monotonic  around  zero  flux,  which  meets  our  expectation.  At  zero 
symmetric  flux,  we  expect  a  switching  of  the  qubit  from  one  classical  state  to  another,  which 
changes  the  magnetic  flux  penetrating  the  readout  SQUID  loop.  However,  the  measured 
curve  does  not  look  like  an  abrupt  switching,  as  we  see  a  slow  change  of  the  flux  instead.  We 
suppose  that  our  antisymmetric  control  line  current  Ia  changes  both  symmetric  and  antisym¬ 
metric  flux,  i.e.  both  the  hight  of  the  potential  barrier  and  the  symmetry  of  the  potential. 
This  we  are  going  to  fix  in  the  next  run  by  changing  the  bonding  scheme  of  the  sample.  Our 
current  bonding  scheme  would  work  with  equal  contact  resistances  at  the  ground  pads  of  the 
antisymmetric  control  line,  which  is  a  rather  naive  assumption.  We  believe  that  the  differ- 
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Figure  6:  Current-voltage  characteristics  of  the  readout  SQUIDs  at  300  rnK. 


ence  in  contact  resistances  causes  unequal  current  distribution  between  the  branches  of  the 
antisymmetric  control  line.  We  are  going  to  change  the  bonding  scheme  in  order  to  have  the 
only  ground  pad  in  the  line,  so  that  the  distribution  of  the  current  in  the  line  will  not  depend 
on  contact  resistances.  Another  possible  reason  for  the  smeared  response  is  the  temperature 
of  this  measurement  of  about  300  mK,  which  is  rather  high. 

1.2.3  New  design  submitted  to  MIT  Lincoln  Lab 

Upon  consulting  with  Dr.  William  Oliver  we  have  slightly  changed  sizes  of  Josephson 
Junctions  in  the  readout  DC  SQUID  to  minimize  the  risk  of  damaging  DC  SQUID  junctions 
during  reactive  ion  etching  fabrication  steps.  We  also  shunted  DC  SQUID  junctions  with 
capacitors  to  reduce  quantum  fluctuations  and  phase  diffusion  in  the  junctions,  thus  the 
shunting  capacitors  should  increase  resolution  of  our  readout  circuit.  The  new  design  was 
submitted  to  MIT  Lincoln  Laboratory  in  June  2008.  We  expect  to  receive  new  chips  from 
MIT  Lincoln  Laboratory  in  November  2008. 

1.3  A1/A10X/A1  junction  based  qubits 
1.3.1  Fabrication  in  Erlangen 

While  waiting  for  chips  ordered  at  MIT  Lincoln  Lab  we  tried  to  realize  the  design  de¬ 
scribed  above  in  Erlangen.  Using  on-chip  biasing  resistors  required  switching  to  a  multilayer 
technological  process.  We  also  had  to  establish  a  process  which  would  allow  us  to  fabricate 
Josephson  junctions  of  size  down  to  0.1  x  0.3  jum2.  It  turned  very  difficult  to  fabricate  such 
junctions  reproducibly  using  our  e-beam  lithography  system  in  Erlangen.  For  this  reason,  we 
obtained  access  to  a  high-end  e-beam  lithography  system  at  the  EU-funded  facility  in  the  de¬ 
partment  of  Microtechnology  and  Nanoscience  (MC2)  at  Chalmers  University  of  Technology 
in  Gothenburg,  Sweden.  This  system  based  on  JEOL  JBX-9300FS  e-beam  machine  is  capable 
of  writing  a  whole  3  inch  wafer  at  once.  Our  approach  was  to  make  all  the  lithography  steps 
using  a  full  wafer  at  MC2,  dice  the  wafer  into  chips  5x5  mm2  in  size,  and  finally  fabricate 
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Figure  7:  Dependence  of  the  critical  current  of  the  readout  SQUID  on  the  chip  B4  from  the 
wafer  DSM4_07_11_4  on  symmetric  control  flux. 

qubits  together  with  readout  circuits  chip  by  chip  by  using  our  aluminum  evaporation  setup 
in  Erlangen.  In  this  procedure,  we  could  still  vary  junction  sizes  and  critical  current  density 
from  chip  to  chip  by  controlling  the  evaporation  angle  and  oxidation  pressure  and  time. 

To  increase  yield  of  working  structures  every  chip  designed  for  MC2  facility  contained 
eight  identical  structures.  One  wafer  contained  100  chips. 

We  divided  the  technological  process  into  four  steps: 

1.  making  alignment  marks; 

2.  depositing  of  palladium  resistors; 

3.  definition  of  contact  pads  and  control  lines; 

4.  fabrication  of  qubits  and  readout  circuits. 

At  every  technological  step  we  used  lift-off  process.  E-beam  lithography  was  performed 
at  acceleration  voltage  of  100 kV  in  order  to  minimize  the  proximity  effect.  For  the  first 
three  lithography  steps,  we  used  combinations  of  LOR  A  and  UV5  resists  because  of  a  small 
exposure  dose  required  for  these  combinations.  Palladium  resistors  were  deposited  using 
electron  beam  evaporation  technique  at  MC2.  To  achieve  better  adhesion,  a  thin  layer  of 
titanium  was  evaporated  before  palladium.  Contact  pads  and  control  line  were  fabricated 
using  thermal  evaporation  of  aluminum  at  MC2.  To  have  a  better  control  of  the  undercut 
during  the  last  lithography  step  we  used  high  resolution  electron  beam  resist  ZEP  520A  as 
top  layer  of  the  resist  sandwich.  The  advantage  of  this  resist  is  that  its  developer  o-Xylene 
does  not  influence  MMA-MAA  resist  that  we  used  as  bottom  layer.  Finally,  we  developed 
MMA-MAA  resist  in  a  mixture  of  isopropyl  alcohol  and  deionised  water.  Unfortunately, 
during  MMA-MAA  resist  prebaking  there  was  an  insulating  polymer  layer  forming  on  top  of 
the  palladium  layer.  We  have  found  that  this  layer  can  be  destroyed  by  passing  a  dc  current  of 
about  0.5  mA  through  the  contact  pads.  However,  this  procedure  occasionally  also  destroyed 
Josephson  junctions  that  we  used  in  the  readout  circuit. 
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Figure  8:  Dependence  of  the  critical  current  of  the  readout  SQUID  on  the  chip  B4  from  the 
wafer  DSM4_07_11_4  on  antisymmetric  control  flux. 


The  SEM  micrographs  of  the  qubit  and  the  readout  circuit  that  we  fabricated  at  MC2 
and  completed  by  aluminum  evaporation  in  Erlangen  are  shown  in  Fig.  [9j  One  can  see  that 
junctions  have  nearly  rectangular  shape.  That  makes  evaluating  of  junctions  areas  easier  and 
more  precise,  thus  helping  to  calculate  qubit  parameters  more  accurately. 


1.3.2  Measurements 


We  tested  structures  fabricated  at  MC2  at  300  mK.  A  current- volt  age  characteristic  of 
a  readout  SQUID  is  presented  in  Fig.  10  One  can  readily  see  a  step  corresponding  to  four 


times  the  energy  gap  of  aluminum  that  is  a  feature  of  three  junction  SQUID. 

A  dependence  of  the  critical  current  of  the  readout  SQUID  on  magnetic  field  is  presented 
in  Fig.  11  From  this  plot  one  can  see  that  the  SQUID  has  non-zero  sensitivity  at  zero 


magnetic  field,  which  in  fact  expected  for  the  used  three  junction  SQUID  design. 

We  also  tested  the  newly  fabricated  structures  in  our  dilution  refrigerator  at  about  25  mK. 
Unfortunately,  we  found  that  the  electromagnetic  noise  was  too  high  to  get  any  reasonable 
results.  Besides  that,  these  structures  turned  out  to  be  extremely  sensitive  to  electro-static 
discharge  (ESD)  and  several  of  them  that  we  cooled  down  we  destroyed.  The  most  critical  step 
for  these  samples  was  removal  of  polymer  layer  between  palladium  resistors  and  aluminum 
electrodes. 

In  order  too  minimize  electromagnetic  interference  and  noise,  we  have  implemented  a  new 
wiring  scheme  in  the  dilution  refrigerator.  We  deliver  currents  to  the  sample  using  single 
ended  coaxial  cables;  the  voltage  drops  will  is  measured  using  shielded  twisted  pairs.  Current 
and  voltage  lines  are  spatially  separated.  This  should  reduce  a  crosstalk  between  voltage  and 
current  lines.  We  also  use  current  dividers  with  ratio  between  100  and  1000  as  well  as  passive 
low-pass  filters  with  cutoff  frequency  around  100  kHz. 
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Figure  9:  SEM  images  of  a  sample  fabricated  at  MC2:  a)  the  two-cell  qubit  with  the  readout 
SQUID;  b)  two  larger  junctions  in  the  readout  SQUID;  c)  inner  qubit  junction;  d)  outer  qubit 
junction. 
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Figure  10:  Current-voltage  characteristic  of  the  readout  SQUID  of  MC2-made  sample. 
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Figure  11:  Dependence  of  the  critical  current  of  the  readout  SQUID  of  MC2-made  sample  on 
magnetic  field. 
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2  High-contrast  readout  of  superconducting  qubits 


2.1  Motivation 

In  the  recent  years,  great  progress  has  been  achieved  in  experiments  with  superconducting 
circuits  operated  in  the  quantum  coherent  regime.  At  millikelvin  temperatures,  these  cir¬ 
cuits  display  quantized  energy  levels  of  which  the  lowest  two  are  used  to  construct  quantum 
bits  [U  [2].  Among  various  qubit  implementations,  Josephson  junction  based  qubits  promise 
an  outstanding  scalability  towards  multi-qubit  circuits  and  are  conveniently  controlled  by  dc- 
and  microwave  frequency  currents.  The  remaining  major  limitation  of  the  superconducting 
approach  to  quantum  computing  is  a  relatively  short  coherence  time  of  the  qubit  states. 

An  obvious  source  of  decoherence  for  superconducting  qubits  is  the  coupling  to  a  detec¬ 
tor  circuit  necessary  for  the  qubit  readout.  A  strongly  coupled  detector  maximizes  readout 
contrast,  while  weaker  coupling  means  a  better  isolation  of  the  qubit  thus  facilitating  longer 
coherence  times.  A  weaker  coupled  detector  disturbs  the  qubit  less,  traded  for  a  loss  of 
readout  fidelity. 

In  this  section,  we  describe  the  development  of  a  data  evaluation  procedure  which  allows 
to  maximize  the  qubit  readout  contrast  also  for  a  weakly  coupled  detector.  We  use  this 
technique  in  experiments  where  we  measure  high-fidelity  Rabi  oscillations  in  Josephson  phase 
qubits.  To  demonstrate  its  usefulness,  we  operate  the  phase  qubit  in  an  elevated  temperature 
regime  where  the  readout  contrast  is  degraded  by  thermal  fluctuations.  We  show  that  by 
using  our  new  data  evaluation  method  the  full  100%  readout  contrast  can  be  recovered. 

2.2  The  Josephson  phase  qubit 

The  logical  states  in  experiments  on  Josephson  phase  qubits  [3]  and  flux  qubits  [3]  are 
distinguished  by  the  magnetic  flux  which  is  generated  by  a  persistent  current  in  a  super¬ 
conducting  loop.  In  all  contemporary  designs  of  these  types  of  qubits,  the  qubit  states  are 
determined  by  a  measurement  of  the  switching  current  of  a  dc-SQUID  used  as  a  flux  detector. 
The  probabilities  of  the  qubit  states  are  thus  obtained  by  a  measurement  of  the  statistical 
distribution  of  the  SQUID  switching  current. 

The  Josephson  phase  qubit  0  consists  of  a  superconducting  loop  interrupted  by  a  Joseph¬ 
son  junction.  The  complete  qubit  circuit  includes  a  superconducting  on-chip  coil  to  provide 
for  flux  bias  and  an  inductively  coupled  dc  SQUID  which  is  used  as  the  detector  for  qubit 
readout.  The  phase  difference  tp  of  the  superconducting  order  parameter  across  the  junction 
depends  on  the  magnetic  flux  threading  the  loop.  The  qubit  loop  inductance  L  is  chosen  such 
that  the  potential  has  the  form  of  a  double  well.  By  changing  the  external  magnetic  flux,  one 
potential  well  is  made  shallower  while  the  other  well  becomes  deeper. 

The  anharmonicity  of  the  shallow  well  allows  to  construct  a  quantum  bit  by  using  the 
two  levels  of  lowest  energy  as  the  logical  states  |0)  and  |1).  Initialization  of  the  qubit  is 
accomplished  by  setting  the  external  magnetic  flux  to  zero,  where  only  one  potential  well 
exists,  and  then  ramping  up  the  flux  adiabatically  to  the  working  point  at  which  only  a 
few  energy  levels  remain  in  the  well.  Logical  qubit  operations  are  accomplished  by  resonant 
photon  absorption  of  an  applied  microwave  field  which  induces  coherent  (Rabi-)  oscillation 
between  the  energy  levels  [3] . 

The  phase  eigenstates  differ  by  the  rate  of  quantum  tunneling  through  the  potential  barrier 
which  separates  the  shallow  and  deep  wells.  Measuring  of  the  qubit  is  accomplished  by  a 
nanosecond-long  pulse  of  magnetic  flux  which  adiabatically  reduces  the  depth  of  the  shallow 
well.  Its  amplitude  is  adjusted  such  that  tunneling  to  the  deep  well  occurs  only  from  the 
excited  qubit  state.  The  readout  pulse  therefore  maps  the  logical  qubit  states  |0)  and  1 1)  to 
the  states  | L)  and  \R)  localized  in  separate  potential  wells.  Since  neighboring  potential  wells 
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correspond  to  states  of  opposite  circulation  direction  of  the  loop  current,  a  transition  from 
the  shallow  to  the  deep  well  results  in  a  change  of  magnetic  flux  [Sj . 

2.3  Switching-current  readout  using  a  dc-SQUID 

The  magnetic  flux  in  the  qubit  loop  is  detected  using  an  inductively  coupled  dc-SQUID 
by  determining  its  switching  current  from  the  superconducting  to  resistive  state.  The  dc- 
SQUID  switching  current  is  measured  by  increasing  its  bias  current  until  a  voltage  appears 
across  the  SQUID  [6].  By  repeating  the  qubit  preparation  and  readout  procedure  several 
thousand  times,  a  statistical  switching  current  distribution  is  acquired.  The  mean  switching 
current  depends  mostly  on  the  SQUID  critical  current,  which  is  a  function  periodic  in  <ho- 
Its  standard  deviation  a  originates  in  the  probabilistic  nature  of  the  switching  process,  and 
is  affected  by  quantum,  thermal  and  bias  current  fluctuations  J7J. 

Distinguishing  the  qubit  states  | L)  and  \R)  by  a  single  shot  switching  current  measurement 
is  only  possible  if  the  histograms  corresponding  to  the  two  states  do  not  overlap.  This  requires 
that  the  difference  between  the  mean  switching  currents  is  about  four  times  larger  than  the 
standard  deviation.  The  change  of  the  qubit  persistent  current  A Ip  ~  Ic  due  to  an  inter-well 
transition  generates  a  flux  change  A<f>sq  =  MAIp  in  the  readout  SQUID.  The  resulting  shift 
of  the  mean  switching  current  depends  on  the  dc-SQUID  bias  flux  and  its  critical  current 
modulation  depth.  The  latter  is  reduced  [8j  to  below  50  %  if  the  dc-SQUID  inductance  Lsq 
and  critical  current  /CjSq  are  large  so  that  LsqIc  ,sq/®0  >  1  • 

To  reach  the  single  shot  readout  limit,  a  typical  minimal  flux  signal  from  the  qubit  should 
be  about  10  m<I>o,  requiring  a  mutual  inductance  of  M  ~  20  pH  at  a  qubit  junction  critical 
current  of  Ic  =  1/xA.  The  mutual  inductance  mediates  dissipation  in  the  dc-SQUID  and 
generates  back-action  to  the  qubit.  While  the  back-action  scales  as  M2/LLsq,  dissipation 
induced  in  the  qubit  is  proportional  to  the  inverse  effective  dc-SQUID  damping  impedance 
Zgq1,  which  is  transformed  to  the  qubit  via  Z~1  =  [M / L)2 .  A  quadratic  benefit  is  hence 

expected  from  a  reduced  mutual  inductance. 

In  practical  qubit  measurements,  statistics  about  the  state  population  are  always  acquired 
because  of  the  probabilistic  nature  of  quantum  mechanics.  A  single  shot  measurement  is 
therefore  not  mandatory  as  long  as  it  remains  possible  to  reconstruct  the  individual  qubit 
state  contributions  from  the  statistical  data  set  with  high  enough  fidelity. 

2.4  Data  evaluation  procedure 

Our  new  method  of  reconstructing  the  qubit  state  population  from  overlapping  SQUID 
switching  current  histograms  is  the  following.  The  qubit  is  first  prepared  in  state  |  L)  and 
the  switching  current  probability  histogram  P(I)  generated  by  this  state  is  recorded.  This 
histogram  is  fitted  to  an  arbitrary  function  //,(/).  For  simplicity,  we  are  using  a  10  degree 
polynomial  truncated  between  minimal  and  maximal  measured  switching  currents  as  the 
fitting  function  /l(/).  An  alternative  fitting  to  a  physically  correct  distribution,  which  is  an 
integral  function  of  the  exponential  escape  rate,  is  also  possible  but  requires  more  computing 
time  and  an  adaption  of  the  model  depending  on  the  sample.  The  same  procedure  is  repeated 
for  a  histogram  obtained  from  measurements  of  a  qubit  prepared  purely  in  the  state  |ii),  thus 
obtaining  the  fitting  function  fR(I)  ■  All  further  overlapping  histograms  are  then  fitted  to  the 
function 

P(I)  =  P(|l»  fR(I)  +  [1  -  P( |1»]  fL(I),  (3) 

hereby  obtaining  P(|l)),  the  population  of  the  excited  qubit  state  1 1) ,  as  the  only  fitting 
parameter. 
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Figure  12:  (a)  Switching  current  distributions  of  a  dc-SQUID  broadened  by  thermal  fluctua¬ 
tions  at  T  =  275  mK.  Left  and  right  histograms  correspond  to  the  qubit  being  in  state  | L) 
and  |  R) ,  respectively.  The  intermediate  histogram  corresponds  to  equal  probability  of  the  two 
states.  Dashed  lines  are  polynomial  fits,  (b)  The  probability  to  measure  the  state  | R)  versus 
externally  applied  flux.  The  dotted  curve  is  calculated  using  the  standard  method,  resulting 
in  a  maximum  contrast  of  70%.  The  solid  curve  has  been  obtained  from  weighted  fits  to  the 
histograms,  which  reconstruct  the  two  states  with  100%  contrast. 


2.5  Experimental  results 

Figure  [l2|a)  shows  histograms  of  the  phase  qubit  prepared  in  states  \L),  |R)  and  in  a 
superposition  state  together  with  the  fitting  functions  /l(7)  and  /l(7)  drawn  as  dashed  lines. 
The  standard  single  shot  method  of  obtaining  the  qubit  population  is  by  setting  a  threshold 
value  /thresh  of  switching  current  and  counting  the  number  of  counts  above  this  threshold 
N(ISW  >  /thresh)  relative  to  the  total  number  of  events  N, 


m))  = 


N(ISW  >  /thresh) 

N  ' 


(4) 


A  loss  of  readout  contrast  here  occurs  because  of  switching  current  events  lying  in  the  region 
of  overlapping  histograms,  which  can  not  be  attributed  to  one  or  other  state.  This  reduced 
contrast  is  shown  in  Fig.  [X2))b)  by  the  dashed  curve  calculated  according  to  Eq.  Q.  The 
alternative  procedure  of  extracting  the  excited  qubit  state  population  P(|l))  by  using  Eq.  ©> 
is  shown  in  Fig.  [I2|)b)  by  solid  line.  One  can  see  that  it  restores  the  full  100%  contrast. 

To  illustrate  the  benefit  of  our  readout  contrast  enhancement  procedure,  in  Fig.  16  we  show 
a  Rabi  oscillation  measurement  with  a  phase  qubit  [9],  provided  to  us  by  the  UCSB  group. 
Here  we  evaluate  the  data  with  both  the  standard  threshold  and  histogram  fitting  methods. 
While  at  the  lower  temperature  the  histograms  are  narrow  enough  to  nearly  reach  the  single 
shot  readout,  at  the  higher  temperature  thermally  broadened  histograms  are  overlapping. 
The  reduction  of  the  Rabi  oscillation  amplitude  due  to  this  overlap  is  fully  eliminated  by 
using  the  fitting  procedure. 


2.6  Application  in  multi-qubit  circuits 

The  technique  demonstrated  above  appears  promising  for  reading  out  multi-qubit  circuits. 
The  number  of  necessary  readout  SQUIDs  can  be  reduced  by  weakly  coupling  several  qubits 
to  a  single  dc-SQUID  and  employing  the  fitting  procedure  to  reconstruct  the  individual  qubit 
flux  states.  After  completion  of  a  quantum  algorithm,  each  qubit  generates  a  flux  signal  in 
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Figure  13:  Rabi  oscillation  in  the  excited  qubit  state,  measured  at  temperatures  of  15  mK 
(top)  and  275  mK  (bottom).  Blue  data  points  (fitted  by  dashed  line)  are  obtained  by  the  stan¬ 
dard  evaluation  method.  A  loss  of  oscillation  visibility  due  to  thermal  histogram  broadening 
is  avoided  by  the  weighted  fit  evaluation,  resulting  in  the  red  data  points  (solid  line). 


the  dc-SQUID  of 


A$sq  =  ^MiA/Cii, 

i 


(5) 


where  A/Cjj  is  the  signal  from  qubit  i  and  Af*  is  the  mutual  inductance  between  qubit  i  and 
the  dc-SQUID.  Our  estimate  shows  that  for  two  qubits  the  ratio  M2 /Mi  must  be  a  factor 
of  two  to  be  able  to  distinguish  the  four  possible  states.  Our  method  allows  to  reduce  the 
coupling  strength  to  the  dc-SQUID  for  all  qubits,  thus  avoiding  decoherence. 


3  Josephson  phase  qubits 

A  vital  factor  in  research  progress  on  macroscopic  quantum  coherent  circuits  is  to  obtain 
an  understanding  of  the  nature  and  origin  of  decoherence.  Obvious  sources  of  decoherence 
that  affect  superconducting  qubits  are  bias  current  and  magnetic  field  fluctuations,  which 
can  be  controlled  by  a  proper  experimental  setup.  Recently,  the  importance  of  the  material 
properties  of  the  Josephson  junction  barrier  was  realized  by  the  UCSB  group  led  by  John 
Martinis,  resulting  in  an  immediate  prolongation  of  the  coherence  time  by  a  factor  of  about 
20.  However,  until  now  this  group  did  no  systematic  experimental  study  of  the  influence 
of  temperature  on  superconducting  qubits.  Common  notion  shared  by  the  superconducting 
qubit  community  is  that  good  qubit  performance  can  only  be  obtained  at  temperatures  below 
100  mK,  whose  generation  requires  expensive  and  bulky  dilution  refrigerators. 

In  contrast  to  the  above  mentioned  common  believe,  our  latest  experimental  results  re¬ 
ported  below  show  that  operation  of  Josephson  phase  qubits  is  possible  at  temperatures  much 
higher  than  100  mK.  We  found  that  the  highest  temperature  where  coherent  Rabi  and  Ramsey 
oscillation  can  be  observed  is  limited  by  the  qubit  energy  level  splitting  which  is  determined 
by  the  circuit  design.  Beneath  providing  additional  hints  towards  a  better  understanding 
of  decoherence  in  superconducting  circuits,  these  new  findings  have  appealing  implications 
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concerning  the  realization  of  solid-state  quantum  computing.  Commercially  available  cryogen- 
free  He3  cryostats  appear  to  be  sufficient  to  cool  the  qubit  chip,  simplifying  both  research 
and  operation  of  quantum  computers. 

As  the  logical  qubit  states  |0)  and  1 1) ,  a  phase  qubit  uses  the  lowest  two  energy  eigenstates 
of  the  Josephson  phase  in  a  metastable  potential  well  of  a  current-biased  junction.  An  elegant 
way  to  decouple  the  junction  from  the  electromagnetic  environment  is  by  applying  its  bias 
current  through  a  dc-transformer  [5],  hereby  embedding  the  junction  in  a  superconducting 
loop  as  it  is  shown  in  Fig.  [l4|a).  The  resulting  circuit  is  known  as  an  rf-SQUID  and  has  a 
potential  energy 

t/(v)  =  ^(1-COS¥’+2k(»’-2,r|f)  )'  (6) 


where  tp  is  the  phase  difference  across  the  junction,  Ic  is  its  critical  current,  <f>ext  the  externally 
applied  flux  through  the  qubit  loop,  and  4>o  is  the  superconducting  flux  quantum.  The  qubit 
loop  inductance  L  is  chosen  such  that  the  parameter  (3l  =  2-kLIc/^q  ~  4,  resulting  in  a 
double  well  potential  plotted  in  Fig.  14  'b).  By  changing  the  flux  4>ext,  one  of  the  wells  can 
be  made  shallow  enough  such  that  it  contains  only  a  small  number  of  energy  levels.  The 
excited  state  |1)  is  populated  by  resonant  photon  absorption  of  an  applied  microwave  field 
which  induces  coherent  Rabi  oscillation  between  the  energy  levels.  Reading  out  the  state 
population  is  done  by  using  a  nanosecond  long  but  still  adiabatic  dc-pulse  of  magnetic  flux, 
which  reduces  the  height  of  the  barrier  separating  the  wells.  The  pulse  amplitude  is  adjusted 
such  that  an  inter-well  transition  occurs  only  from  the  excited  state.  Since  each  potential 
minimum  corresponds  to  a  certain  circulation  direction  of  the  loop  current,  the  transition 
from  the  shallow  to  the  deep  well  results  in  a  change  of  magnetic  flux,  which  is  afterwards 
read  out  by  recording  the  switching-current  distribution  of  an  inductively  coupled  dc-SQUID. 


3.1  Samples 

We  studied  the  influence  of  temperature  on  qubit  performance  using  two  different  samples. 
Sample  #1  was  fabricated  according  to  our  designs  at  VTT  foundry  [TO]  using  a  standard 
lithographic  Nb/A10x/Nb-trilayer  processes,  featuring  Josephson  junctions  of  7  ^m2  mini¬ 
mum  area  and  30  A/cm2  current  density.  Sample  #  2  has  been  fabricated  by  John  Martinis 
group  at  UCSB  in  a  custom  process,  using  aluminum  as  the  electrode  material  and  a  small 
qubit  junction  of  size  approximately  1  /xm2.  In  contrast  to  the  standard  VTT  samples  fea¬ 
turing  Si02  insulation,  in  the  UCSB  sample  SiNx  is  used  as  the  dielectric. 


(a) 


(b) 


Figure  14:  (a)  Schematic  of  the  phase  qubit  circuit.  The  qubit  junction  is  embedded  in  a 
superconducting  loop  which  is  coupled  inductively  to  the  flux-biasing  coil  and  readout  dc- 
SQUID.  (b)  Sketch  of  the  qubit  potential  U(ip).  On  the  right  side,  a  zoom  into  the  shallow 
left  potential  well  indicates  the  wave  functions  describing  the  two  qubit  states  |0)  and  |1). 
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Figure  15:  The  widths  5En  of  the  energy  levels  n  are  indicated  by  grey  bars  and  the  energy 
level  separations  are  shown  by  black  arrows. 


The  UCSB  group  has  demonstrated  recently  m  that  minimizing  the  junction  area  and 
thus  the  number  of  two-level  fluctuators  in  the  junction  tunnel  barrier  significantly  improves 
qubit  performance.  Tuning  the  junction  capacitance  in  the  appropriate  range  is  achieved 
by  using  a  silicon  nitride  capacitive  element  separated  on  chip  from  the  Josephson  tunnel 
junction.  SiNx  does  not  suffer  from  two-level  fluctuators  and  thus  has  much  lower  dielectric 
loss  compared  to  the  oxide  dielectrics  Si02-  Accordingly,  measurements  of  the  UCSB  sample 
#2  in  our  setup  showed  coherence  times  which  were  by  more  than  one  order  of  magnitude 
longer  than  in  the  Nb  SiCU  based  sample  #1. 


3.2  Spectroscopic  quality  factor 

To  measure  the  life  time  T\  of  the  excited  state  we  apply  a  resonant  microwave  7r-pulse 
which  prepares  the  qubit  in  the  excited  state.  Then,  after  a  variable  delay  time,  we  measure 
the  remaining  population  in  the  excited  state.  We  observed  exponential  decay  at  half-life 
times  of  about  4  ns  for  sample  #1  and  about  90  ns  for  sample  #2. 

Quantum  mechanically,  the  life  time  Tj  is  related  to  an  uncertainty  in  the  energy  of  the 
excited  state,  giving  rise  to  a  broadening  5E  of  the  level  at  energy  E  as  it  is  sketched  in 
Fig.  15  Therefore,  in  consistence  with  the  measured  Tj  relaxation  times,  we  observed  that 


the  spectroscopic  resonance  peak  for  photon  activation  to  the  excited  state  has  a  width  in 
frequency  6lu  of  about  27r-  10  MHz  for  sample  #2  and  about  27T-400  MHz  for  sample  #1. 
Corresponding  Q  factors  calculated  by 


=  ^  =  T 

^  6i0  5E 


(7) 


were  about  30  at  a  transition  frequency  cuoi  =  27t  •  16.5  GHz  for  the  sample  #1  and  about 
600  at  a  frequency  woi  =  27r  •  8  GHz  for  sample  #2.  The  anharmonicity  of  the  qubit  potential 
has  the  consequence  that  the  separation  between  adjacent  levels  decreases  with  increasing 
energy.  An  external  microwave  will  hence  be  resonant  only  to  one  transition,  provided  that 
the  individual  levels  are  not  too  broad.  From  Fig.  15  it  can  be  seen  that  this  corresponds  to 
the  condition 

1"~  '  (8) 


+  -(5 Ei  +  6E2)  <C  huj oi, 


where  SEn  is  the  full  width  of  level  n  arising  from  its  finite  lifetime  T”1  and  iomn  is  the 
transition  frequency  from  level  m  to  level  n. 

It  is  important  to  note  that  the  qubit  potential  well,  as  its  depth  is  increased,  becomes 
more  harmonic  and  the  difference  in  the  transition  frequencies  between  neighboring  levels 
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gets  smaller.  Taking  into  account  the  broad  resonance  line  width  of  Nb-based  qubit  #1,  we 
can  therefore  not  regard  this  qubit  as  a  two-level  system,  but  should  rather  expect  multilevel 
dynamics  m-  In  contrast,  sample  #2  is  operated  by  applying  microwaves  as  a  true  two-level 
system  due  to  its  high  Q  value.  Deeper  potential  well  thus  remains  anharmonic  enough  not 
to  result  in  populating  of  higher  energy  levels  by  microwave  applied  at  Uq\. 


3.3  Coherent  oscillation  at  high  temperatures 


We  observed  Rabi  oscillation  in  the  population  of  the  excited  state  by  varying  the  duration 
of  a  microwave  pulse  which  frequency  is  resonant  to  the  |0)  to  |1)  transition.  By  fitting  the 
data  to  an  exponentially  damped  sine  function,  we  extract  Rabi  amplitude,  frequency  and 
decay  time.  The  above  reported  Nb-based  phase  qubits  show  rather  modest  decoherence 
times,  where  samples  from  VTT  perform  very  similarly  to  our  previously  measured  Hypres 
samples.  In  these  Nb-based  qubits  Ti  as  well  as  Rabi  decay  time  Tr  were  typically  3  to  7  ns, 
which  is  significantly  shorter  than  the  data  for  Al-based  circuits. 

In  order  to  observe  Rabi  oscillation  at  higher  temperatures,  we  have  to  avoid  thermal 
activation  out  of  the  shallow  potential  well.  The  activation  rate  increases  as  the  potential 
barrier  height  becomes  comparable  to  the  thermal  energy  ksT.  By  reducing  the  field  bias 
and  operating  in  a  deeper  potential  well,  we  were  able  to  observe  oscillations  at  temperatures 
up  to  0.9  K  in  the  Nb-based  qubit  #1.  At  a  temperature  of  0.8  K,  the  amplitude  was  reduced 


by  a  factor  of  one  half,  while  the  oscillation  decay  time  dropped  only  by  80%,  see  Fig.  17 


We  found  that,  above  a  certain  temperature,  it  becomes  necessary  to  reduce  the  amplitude 
of  the  readout  pulse  in  order  to  preserve  good  contrast  of  oscillations.  This  is  due  to  the  fact 
that,  during  the  readout  pulse,  activation  from  the  excited  qubit  state  to  the  deep  well  becomes 
possible  not  only  by  quantum  tunnelling  but  also  by  thermal  activation  over  the  barrier. 
Moreover,  at  higher  temperatures  the  width  of  the  measured  switching-current  histogram 
of  the  readout  dc-SQUID  increases.  We  were  able  to  avoid  a  loss  of  contrast  between  the 
two  qubit  states  due  to  broadened  histograms  by  a  fitting  procedure  extracting  the  additive 
weighted  contributions  of  the  two  states. 

Our  results  for  the  high-Q  sample  #2  are  summarized  in  Fig.  p~8|  We  plotted  here  the 
observed  dependence  of  the  Rabi  oscillation  amplitude  and  decay  on  temperature.  Remark¬ 
ably,  we  notice  only  a  very  small  decrease  of  decoherence  times  at  low  temperatures.  The 
oscillations  disappear  completely  close  to  the  temperature  T  =  fiu/kB,  which  is  355  mK  at 
ui  =  27 r-  7.4  GHz  and  384  mK  at  oj  =  27t-  8  GHz. 

To  measure  the  influence  of  temperature  on  the  dephasing  time,  in  sample  #2  we  per¬ 
formed  a  Ramsey- type  experiment,  in  which  two  7r/2-pulses  separated  by  a  variable  duration 
are  applied  to  the  qubit.  As  expected,  the  frequency  of  the  observed  Ramsey  fringes  was 
equal  to  the  detuning  of  the  microwave  pulses.  We  took  a  set  of  data  at  335  MHz  below  the 


exact  resonance  at  8  GHz.  Figure  19  shows  the  Ramsey  oscillation  amplitude  and  extracted 
T2  time  versus  temperature.  These  oscillations  decay  similarly  to  the  Rabi  data. 


4  Summary  and  conclusion 

We  redesigned  our  structures  to  realize  the  newly  proposed  “silent”  two-cell  qubits.  At 
the  operating  point,  the  “silent”  qubit  should  be  insensitive  to  fluctuations  caused  by  the 
flux  control  circuits.  We  optimized  our  design  for  different  critical  current  densities  and 
submitted  it  to  MIT  Lincoln  Laboratory.  We  have  observed  a  switching  of  the  two-cell 
flux  qubit  with  a  passive  7r-shifter  between  two  classical  states  during  a  characterisation 
experiment  of  chips  fabricated  at  MIT  Lincoln  Laboratory.  The  qubit  switches  around  zero 
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Figure  16:  Rabi  oscillations  observed  in  Nb-SiOx  based  qubit  sample  #1  (a)  and  Al-SiN^ 
based  qubit  sample  #2  (b),  at  the  indicated  temperatures.  Solid  lines  are  a  fits  to  an  expo¬ 
nential  decaying  sine  functions  from  which  Rabi  amplitude  and  decay  time  are  extracted. 
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Figure  17:  Rabi  oscillation  amplitude  (top)  and  decay  time  (bottom)  versus  temperature, 
measured  using  sample  #1.  Each  data  set  correspond  to  an  adjusted  magnetic  field  resulting 
in  resonance  at  the  frequency  indicated  in  the  legend. 
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Figure  18:  Amplitude  (right  axis)  and  decay  time  (left  axis)  of  Rabi  oscillation  observed  in 
sample  #2,  plotted  versus  temperature.  Each  set  of  data  points  corresponds  to  an  adjusted 
magnetic  field  bias  resulting  in  an  energy  level  spacing  causing  resonance  at  the  microwave 
frequency  indicated  in  the  legend.  Dashed  lines  are  guides  to  the  eye. 
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Figure  19:  Ramsey  oscillation  amplitude  (right  axis)  and  decay  time  (left  axis)  vs.  tempera¬ 
ture  (sample  #  2).  Dashed  lines  are  guides  to  the  eye. 


flux  for  both  symmetric  and  antisymmetric  flux  components,  which  meets  the  expectation  of 
our  experiment.  This  behaviour  minimises  the  coupling  of  the  qubit  to  the  environment. 

We  developed  a  new  multilayer  technological  process  consisting  of  fabrication  of  the  struc¬ 
tures  partially  at  MC2  foundry  at  Chalmers  University  and  finishing  the  fabrication  in  Erlan¬ 
gen.  Using  this  process,  we  fabricated  aluminum  junctions  as  small  as  0.09  x  0.24  /mi2,  which 
meets  the  goal  set  for  our  experiment.  We  have  revised  noise  issues  and  wiring  in  our  dilution 
refrigerator.  These  measures  helped  us  to  reduce  interference  and  channels  cross-talk  during 
measurements. 

In  the  second  part  of  the  report  we  presented  a  new  method  of  qubit  data  evaluation  which 
makes  it  possible  to  enhance  the  readout  contrast  in  experiments  with  superconducting  qubits. 
This  technique  employs  a  weighted  fitting  of  switching-current  histograms  obtained  from  the 
dc-SQUID  detector,  which  is  used  for  qubit  readout.  We  demonstrated  reconstructing  of 
the  full  contrast  of  Rabi  oscillation  using  a  phase  qubit  in  a  regime  where  the  oscillation 
visibility  was  reduced  due  to  thermal  fluctuations  in  the  environment.  Similarly,  a  large 
readout  contrast  can  be  obtained  in  qubit  circuit  designs  where  the  detector  SQUID  is  weekly 
coupled  to  the  qubit.  This  feature  is  of  particular  importance,  because  for  achieving  long 
coherence  times  the  qubits  should  be  well  isolated  from  the  environment. 

The  demonstrated  data  evaluation  procedure  will  be  beneficial  for  multi-qubit  circuits, 
offering  an  opportunity  to  read  out  the  states  of  several  qubits  by  a  single  dc-SQUID.  We 
expect  that  it  will  be  thus  possible  to  reconstruct  the  individual  qubit  states  with  good  fidelity. 

In  the  last  part  of  the  report  we  presented  Rabi  oscillations  in  Josephson  phase  qubits 
fabricated  using  standard  Nb  process  with  SiCR  on-chip  dielectric.  Operating  in  a  deep 
potential  well,  Rabi  like  oscillation  persisted  at  temperatures  up  to  0.9  K,  where  only  modest 
decrease  in  decay  time  was  observed.  These  results  are  obtained  in  the  multi-level  limit. 
In  contrast,  A1  based  phase  qubit  with  SiNa,  dielectric  in  the  shunting  capacitor  which  we 
obtained  from  UCSB  has  long  coherence  times  and  can  be  treated  as  a  true  two-level  system. 
Decoherence  times  measured  in  the  UCSB  sample  are  nearly  independent  of  temperature  up 
to  the  point  at  which  the  thermal  energy  becomes  comparable  to  the  energy  level  separation. 
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These  results  suggest  that  presently  existing  best  phase  qubits  designed  to  have  a  high  Q  and 
large  energy  level  separation  could  be  easily  operated  at  higher  temperatures. 


References 

[1]  J.  Q.  You  and  F.  Nori,  Phys.  Today  58(11),  42  (2005). 

[2]  G.  Wendin  and  V.  Shumeiko,  in  Handbook  of  Theoretical  and  Computational  Nanotech¬ 
nology ,  edited  by  M.  Rieth  and  W.  Schommers  (ASP,  Los  Angeles,  2006). 

[3]  J.  M.  Martinis,  S.  Nam,  J.  Aumentado,  and  C.  Urbina,  Phys.  Rev.  Lett.  89,  117901 

(2002). 

[4]  I.  Chiorescu,  Y.  Nakamura,  C.  J.  P.  M.  Harmans,  and  J.  E.  Mooij,  Science  299,  1869 
(2003). 

[5]  R.  W.  Simmonds,  K.  M.  Lang,  D.  A.  Hite,  D.  P.  Pappas,  and  J.M.  Martinis,  Phys.  Rev. 
Lett.  93,  077003  (2004). 

[6]  T.  A.  Fulton  and  L.  N.  Dunkleberger,  Phys.  Rev.  B  9  ,  4760  (1974). 

[7]  A.  Wallraff,  A.  Lukashenko,  C.  Coqui,  A.  Kemp,  T.  Duty  and  A.  V.  Ustinov,  Rev.  Sci. 
Inst.  74,  3740  (2003). 

[8]  B.  Chesca,  R.  Kleiner  and  D.  Koelle  in  The  SQUID  Handbook,  Vol.  I:  Fundamentals 
and  technology  of  SQUIDs  and  SQUID  systems.  Eds.  J.  Clarke  and  A.  I.  Braginski, 
Wiley-VCH,  Weinheim,  (2004)  p.  43. 

[9]  J.  Lisenfeld,  A.  Lukashenko,  M.  Ansmann,  J.  M.  Martinis,  and  A.  V.  Ustinov,  Phys. 
Rev.  Lett.  99,  170504  (2007). 

[10]  VTT  Technical  Research  Center  of  Finland. 

[11]  M.  Steffen,  M.  Ansmann,  R.  McDermott,  N.  Katz,  R.  C.  Bialczak,  E.  Lucero,  M.  Neeley, 
E.  M.  Weig,  A.  N.  Cleland,  and  J.  M.  Martinis,  Phys.  Rev.  Lett.  97,  050502  (2006). 

[12]  J.  Claudon,  F.  Balestro,  F.  W.  J.  Hekking,  and  O.  Buisson,  Phys.  Rev.  Lett.  93,  187003 
(2004). 


22 


